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Abstract. A new data presentation is proposed in order to consider anomalous J/'P suppression in Pb + Pb 
collisions at P/A = 158 GeV/c. If the inclusive differential cross section with respect to a centrality variable 
is available, one can plot the yield of J/U' events per Pb-Pb collision as a function of an estimated 
squared impact parameter. Both quantities are raw experimental data and have a clear physical meaning. 
As compared to the usual J/<I' over Drell-Yan ratio, there is a huge gain in statistical accuracy. This 
presentation could be applied advantageously to many processes in the field of nucleus-nucleus collisions 
at various energies. 



1 Introduction 

Very interesting results have been obtained recently by the 
NA50 experiment at CERN concerning J/'F production 
in Pb + Pb collisions at P/A = 158GeV/c |i|,|J3[|,|. In 
most central collisions, the J/'F events are significantly 
suppressed with respect to what is expected from normal 
nuclear absorption as measured in lighter systems [pp0,p| . 
According to theoretical predictions made more than ten 
years ago by Matsui and Satz Q, this anomalous J/F 
suppression could be a sign of the awaited formation of a 
quark-gluon plasma in nucleus-nucleus collisions at very 
high energy. 

The deficiencies of customary data presentations, us- 
ing the ratio between J/'F and Drell-Yan events or the 
differential cross section for J/F production with respect 
to some measured centrality variable, will be stressed first. 
A new data presentation 0] will then be proposed, which 
removes the previous deficiencies with the help of one key 
additional ingredient, the inclusive differential cross sec- 
tion with respect to the centrality variable. For any pro- 
cess "p" , the yield of "p" events per nucleus-nucleus col- 
lision is plotted as a function of an estimated squared im- 
pact parameter. This new presentation will be applied to 
NA50 results from their 1995 data taking, available in a 
thesis Q . Implications of this new presentation will also 
be discussed. Finally, conclusions will be drawn in the 
perspective of RHIC and LHC experiments on nucleus- 
nucleus collisions at very high energy. 
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2 Usual data presentation 

In the usual presentation of NA50 results concerning J/'F 
production in nucleus-nucleus collisions ||l|,p|,pl0] , the ratio 
between J/'F and Drell-Yan events is plotted as a function 
of the transverse energy E^ measured in an electromag- 
netic calorimeter. This centrality variable is aimed pri- 
marily at sorting out all events according to the impact 
parameter of the collision. The Drell-Yan process, sup- 
posed to be insensitive to nuclear matter effects, is in- 
deed a good reference for normalizing J/'F events from 
the physics point of view. Moreover, some systematic ef- 
fects cancel in such a ratio. However, when one sees an 
interesting feature in a ratio, it is not always obvious to 
know whether it is due to the numerator or the denomi- 
nator. More importantly, since the Drell-Yan continuum 
is much less populated than the J/'F peak in the dimuon 
mass spectrum, the statistical uncertainty on the ratio 
comes essentially from the denominator. It is between 5 
and 10 times larger - 5 for central collisions, 10 for pe- 
ripheral ones - than the contribution from the number of 
events in the J/F peak. In other words, one would need 
between 25 and 100 times less running time, all other con- 
ditions staying equal, to get a given relative statistical un- 
certainty on J/tf' production from the number of events 
in the J/'F peak than from the ratio between J/F and 
Drell-Yan events. This is a first deficiency of the usual 
presentation. One would like to use another quantity for 
normalizing J/'F events without losing so much in statis- 
tical accuracy. 

The ratio between J/'F and Drell-Yan events is ob- 
tained from the differential Et distributions of cross sec- 
tion for both classes of events, which are the basic experi- 
mental data one has to start with. These raw experimental 
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data, da^jdE'Y and i?T, where "p" stands for either J/tf' 
or Drell-Yan process, do not have a very direct physical 
meaning. From the increase or the decrease of dcp/dx as a 
function of -Et , one cannot even infer whether the produc- 
tion of "p" events increases or decreases with the centrality 
of the coUision. In particular, the decrease of Ao^JAEt at 
high E'Y for any "p" process simply reflects the fact that, 
for most central collisions, there is a maximum value of -Et 
beyond which there is no more cross section. Et is surely 
increasing with centrality, but at which rate? The answer 
to this question is needed if one wants to go from Et to 
a more direct centrality variable like the impact param- 
eter. For all these reasons, dcp/dE^T and E^t are rather 
difficult to understand and compare directly with simple 
models. This is the second deficiency of the usual presenta- 
tion. One would like to find other experimental quantities, 
not too far from dup/dET and Et, which would have a 
more direct physical meaning, from which one could di- 
rectly say something on the variation of the production of 
"p" events with centrality, and which could be compared 
directly with simple models. 



3 New data presentation 

One key quantity that could be used to remove the above- 
mentioned deficiencies is the inclusive distribution of the 
cross section for the centrality variable, denoted as C here- 
after for more generality. This inclusive distribution will 
be needed in its differential form ddinc/dC for normaliz- 
ing dcTp/dC, and in its integral form (Tinc(C') for getting 
an estimated squared impact parameter {\P')c- 

When we use dcp/dC we mix the probability to get a 
given centrality with the probability to get the "p" pro- 
cess at this centrality. This is precisely why the increase or 
the decrease of dcTp/dC as a function of centrality has no 
straightforward meaning. This distribution dcTp/dC is in 
fact the product of two quantities which themselves have a 
more direct physical meaning than their product. It can be 
written as Fp • do'inc/dC', where the inclusive distribution 
duinc/dC carries the probability that a nucleus-nucleus 
collision occurs at a given value C of the centrality vari- 
able, and Kp is the yield of "p" events per nucleus-nucleus 
collision at this given centrality. This yield Y^^ which, be- 
ing equal to (dtTp/dC)/(dcrinc/dC), is a well-defined phys- 
ical quantity. For copiously produced particles it is simply 
their average multiplicity per nucleus-nucleus collision at 
a given centrality. As it is a ratio, it should be insensitive 
to some systematic uncertainties. Its variation as a func- 
tion of C should accurately reflect whether the production 
of "p" events increases or decreases with the centrality 
of the collision. It should not be subject to any artifi- 
cial decrease for most central collisions. For both J/iZ' and 
Drell-Yan processes, one expects that this yield steadily 
increases towards more central collisions, like the number 
of nucleon-nucleon collisions they originate from, unless 
the J 1'^ is very strongly suppressed. The first step in the 
new data presentation is thus to use the yield Fp of "p" 
events per nucleus-nucleus collision instead of dcTp/dC. 



The idea behind tagging a process with a centrality 
variable in nucleus-nucleus collisions is always to sort out 
events according to the impact parameter. If a centrality 
variable C is assumed to vary monotonically as a func- 
tion of the impact parameter 6 - and centrality variables 
are purposely chosen for that reason -, it is very easy 
to go from C to an estimate of &, or more precisely &^. 
One only has to use the integral inclusive cross section 
CinclC*), from most central collisions to any given value 
of C. From the geometrical dependence of the inclusive 
cross section, dcrinc = 27r-& ■ Ah = -k ■ d(6^), one simply gets 
CinclC) = 7r(6^)e where (6^)e is an estimate of the squared 
impact parameter corresponding to the value of C. There 
is a one-to-one correspondence between C and (6^)c. C 
slices are transformed into (&^)c slices with a width pro- 
portional to the number of counts in the C slices. For this 
reason, which is also related to the fact that d(Tinc/d& goes 
to zero at zero impact parameter, (&^)o seems a better 
variable than the estimated impact parameter 6e- Instead 
of dividing crinc(C') by tt, one could divide it by the geo- 
metrical cross section a^^o and get a quantity proportional 
to (6^)e, but with such a normalization that it varies be- 
tween and 1 from most central to most peripheral col- 
lisions. Another quantity which could be interesting to 
use for plotting results from different systems would be 
(^max ~ (^^)e), where &^ax = o'goo/Tr. It has the advantage 
of being correlated, and not anticorrelated, with the cen- 
trality, and extends to larger and larger values for larger 
and larger systems, with the zero value corresponding al- 
ways to most peripheral collisions. Such a transformation 
from (Tinc(C') to (&^)c has been used in more or less details 
by several experiments in the field of nucleus-nucleus col- 
lisions at various incident energies. Its reliability has been 
discussed thoroughly, and has been checked to be excellent 
within the framework of the intranuclear cascade model at 
energies per nucleon around 1 GeV |0] . Model calculations 
are obviously needed to evaluate the method and to com- 
pare the quality factors of various centrality variables Q , 
which combine the fluctuations of C at any given h and 
the variation rate of C with respect to 6 or \? . For some 
AGS or SPS experiments, even though the data are plot- 
ted as a function of a centrality variable, a scale for the 
impact parameter estimated along the preceding lines is 
indicated in parallel (ll]jlj] . Finally, the second step in the 
new data presentation consists in replacing the measured 
centrality variable C with an estimate of the squared im- 
pact parameter, (fe^)e = crinc(C')/7r. 

After applying both steps one gets the yield Fp of "p" 
events per nucleus-nucleus collision as a function of the 
estimated squared impact parameter. Both quantities are 
raw experimental data and have a clear physical meaning. 
As compared with the usual ratio between J/'P' and Drell- 
Yan events, there is a huge gain in statistical accuracy for 
J/if. 

With this new data presentation one can consider any 
process independently of all others, with the best statistics 
available for each of them. In the same way as the integral 
of dcTp/dC as a function of C is the total cross section for 
the "p" process, the integral of Fp as a function of {6^)e 
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is equal to the total cross section for the "p" process di- 
vided by TT. There is no loss of information in going from 
the centrality variable C to the estimated squared impact 
parameter (5^)e- The limits of the slices used for looking at 
the variation with centrality have simply to be specified for 
both C and (6^)e- The main requirement for this new pre- 
sentation is a good inclusive centrality distribution, with 
high enough statistical accuracy and proper corrections 
for efficiency and empty-target contribution. The yield Yp 
being the ratio of cross sections, part of systematic uncer- 
tainties are removed if inclusive measurements are taken 
simultaneously with the measurements of the "p" process. 
The normalization uncertainty of the inclusive cross sec- 
tion has an effect on the abscissa rather than on the ordi- 
nate, which may be unusual but does not bring about any 
practical problem. 

Such a presentation provides an excellent starting point 
for comparisons between experiment and theory, and also 
between experimental results themselves. Since there is no 
explicit appearance of C in the new presentation, results 
obtained with various centrality variables should be identi- 
cal as long as these variables sample the impact parameter 
the same way. Anyway, such a comparison could help to 
check systematic uncertainties. For comparisons between 
experiment and theory, it is straightforward to compare 
the measured yields as a function of the estimated 5^ with 
the calculated ones as a function of the real h^ . This is 
particularly interesting for a quick comparison with sim- 
ple models. However, in order to take into account the 
fluctuations of any centrality variable as a function of the 
impact parameter, a comparison with better quality would 
result from using the same procedure of 6^ estimation for 
both experiment and theory, even if inclusive cross sec- 
tions do not agree within a high degree of accuracy , or 
from unfolding the experimental results from these fluctu- 
ations. It is also clear that this whole presentation could 
be applied advantageously to other processes than J j'^ 
and Drell-Yan production. 



4 Application 

In a thesis by Bellaiche Q from the NA50 collaboration, 
all necessary pieces of information are available from the 
1995 data taking for applying this new data presentation. 
They have not been published as such. The results pre- 
sented below thus have to be considered with care. They 
are only indicative, and they are to be used simply as an 
illustration of the advantages inherent to the new data pre- 
sentation. Basic experimental data are the differential E^ 
distributions of the cross section for the J/!f' and Drell- 
Yan events (Fig. ^). The key additional ingredient is the 
differential E't distribution of the inclusive cross section, 
by the integration of which one can estimate the squared 
impact parameter (6^)c for each value of fi'x (Fig. ||). Un- 
certainties on dcTinc/d-ET have been neglected in the fol- 
lowing. The inclusive cross section is only available with 
arbitrary units in the thesis. A normalization factor had 
to be introduced to get the (6^)c values in units of fm^. It 
has been adjusted in such a way that the dependence of 
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Fig. 1. Differential transverse-energy distribution of cross sec- 
tion for J 1'^ (top) and Drell-Yan (bottom) production in Pb- 
Pb collisions at PjA = 158 GeV/c § 



(6^)c upon E'Y: agrees with the correlation between the av- 
erage values of i?T and h listed in NA50 publications [QJg] 
for successive E^ slices, as fitted on the basis of a Glauber 
model calculation. These average values are also shown in 
the bottom part of Fig. |[ E'^ values from p have been di- 
vided by 0.74 to take into account the different E'Y scales 
used in the NA50 publications. 

After the first step, i.e. the normalization of the J j'^ 
and Drell-Yan cross sections to the inclusive one, one gets 
(Fig. H) the yields of J/'f' and Drell-Yan events per Pb- 
Pb collision as a function of i?T • Both yields increase with 
i?T, without any artificial decrease at large E^y- Whereas 
this increase is rather steady for Drell-Yan events, there 
is a change of behaviour for J j'^ at about 40 GeV, an E^ 
value beyond which the increase is definitely slower for 
J I'i' than for Drell-Yan events. This is an indication for 
J I'i' suppression in central collisions, relative to Drell-Yan 
events. 

After the second step, i.e. replacing E^ by crinc(-BT)/7r, 
one gets (Fig. g) the yields of J/tf' and Drell-Yan events 
per Pb-Pb collision as a function of (6^)o, the squared 
impact parameter estimated from the E^ inclusive cross 
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Et (GeV) 

Fig. 2. Top: differential transverse-energy distribution of in- 
clusive cross section in Pb-Pb collisions at PjA — 158 GeV/c 
Ig]. Bottom: estimated squared impact parameter versus trans- 
verse energy; the correspondence between Et and b from ||l|J 
is indicated with closed and open circles, respectively 



section. Slices with almost constant width in (6^)e have 
been used. The same remarks could be made as from 
Fig. I after the first step. The Et value of 40 GeV for 
the change of behaviour of the J/'F yield is changed into 
a (6^)e value of 80 fm^. Points with large error bars at 
large Et in Fig. H are all contained in the point at the 
smallest value of (&^)o in Fig. |j. The limits of the yields 
for most central collisions are more easily readable from 
Fig. U than from Fig. g. They could be directly compared 
with the J/S' and Drell-Yan yields in p-p collisions - J/'F 
and Drell-Yan cross sections divided by the total inelastic 
p-p cross section - times the number of nucleon-nucleon 
collisions in most central Pb-Pb collisions from a Glauber 
model calculation. From integration of the yields in Fig. || 
one can get the total cross sections for J/F and Drell-Yan 
production divided by tt. 

A comparison is also made in F ig.H with a model calcu- 
lation a la Blaizot and OUitrault jiy The yield of Drell- 
Yan events per Pb-Pb collision as a function of (6^)c is 
compared, within a scale factor, to the number of nucleon- 



100 120 140 160 180 200 

ET(GeV) 

Fig. 3. Yields of J/^ (top) and Drell-Yan (bottom) events 
per Pb-Pb collision versus transverse energy, calculated from 
NA50 data at P/A = 158 GeV/c ||| 



nucleon collisions calculated in the Glauber model as a 
function of the real 5^, without taking into account the 
fluctuations between the estimated and actual b^. The 
agreement is reasonable. For J/W events this number of 
nucleon-nucleon collisions is multiplied by two correction 
factors for absorption. The first one corresponds to normal 
J/F absorption in nuclear matter with some cross section 
CTabs- The second one is intended for simulating complete 
J/F suppression due to quark-gluon plasma formation. 
It goes down from one to zero when nucleon-nucleon col- 
lisions producing J/F occur in a tube of nuclear matter 
with nucleon density per unit area larger than a critical 
value Pcrit- With (Tabs — 6.0 mb and Pcrit = 2.9 fm^^, the 
model reasonably accounts for the experiment, in particu- 
lar for the clear change of behaviour at an impact param- 
eter of about 9 fm. 

Finally, very accurate results are obtained for J/F pro- 
duction which can be compared easily with simple mod- 
els. The onset of the anomalous J/F suppression can be 
looked at with much better accuracy than on the basis of 
the usual J/F over Drell-Yan ratio. However, we want to 
recall the word of caution from the beginning of this sec- 
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estimated b^ (fm^) 

Fig. 4. Yields of J/iP' (top) and Drell-Yan (bottom) events per 
Pb-Pb collision versus estimated squared impact parameter, 
calculated from NA50 data at PjA = 158 GeV/c |8| (crosses) 
and from the number of nucleon-nucleon collisions within a 
Glauber model (solid lines), with correction factors for J/<Z' 
absorption |1^ (see text for more details) 



tion. Definite conclusions about J/'l' anomalous suppres- 
sion need to be drawn from official data. This is also why 
there was no attempt to calculate error bars for the val- 
ues of CTabs and Pcrit- Morcovcr the Drell-Yan production 
remains an essential result. One has to check its normal 
behaviour within its inherently limited accuracy. 



5 Discussion 

One idea from this new data presentation, the normaliza- 
tion of J/tf' events to the inclusive, or minimum bias, E^ 
distribution, has been used recently by the NA50 collab- 
oration 1^ , making the most of the whole statistics avail- 
able for J/'F production in their 1996 data taking. How- 
ever, this new presentation is not used as such, except the 
first step for Drell-Yan production only. For easy com- 
parison with previously published results, the J/W yield 
is transformed into a "minimum bias" J/'F over Drell- 
Yan ratio, through a division by a model calculation for 



the Drell-Yan yield. In order to stick more closely to the 
raw experimental data, it would be very interesting if the 
new presentation were to be applied as a whole to these 
most recent and also to future NA50 data. One would 
not have to worry anymore because of the different E^ 
scales used in successive presentations. More importantly, 
it would be particularly helpful to compare between one 
another the results obtained with the three centrality vari- 
ables available in this experiment, namely the transverse 
energy Et measured in an electromagnetic calorimeter, 
the zero-degree energy measured in a hadronic calorime- 
ter, and the multiplicity measured in a silicon detector. 
Perhaps it would also be possible to get more accurate 
information on J/W production from older data takings, 
for example in S -I- U collisions. 

Finally, it is interesting to try and quantify the gain 
brought about by the normalization to the inclusive cen- 
trality distribution in the assessment of the anomalous 
J/'F suppression in Pb -I- Pb collisions. It can be done for 
example on NA50 results as shown in Fig. 9 from S. In 
this figure, the "minimum bias" as well as the measured 
J/tf' over Drell-Yan ratios are divided by the normal ab- 
sorption factor and plotted as a function of the mean nu- 
clear path length L (Fig. |5|). The normal absorption ap- 
pears as a horizontal line at a constant value of 1, without 
any information on its uncertainty. In Fig. 0, an uncer- 
tainty band, necessary for a quantitative comparison to 
Pb -I- Pb results, has been added around the straight ref- 
erence line. It has been calculated from the same p nucleus 
and S -I- U data as used in Q , with correct error bars as 
compared to previous NA50 publications (see note added 
in proof to P|), and taking into account the correlation 
between the normalization and the slope of the exponen- 
tial fit. An uncertainty band had already been shown in 
[ [l4[ but it had been calculated with the old error bars for 
S -|- U data and without taking into account the correla- 
tion between the normalization and the slope. By chance 
this uncertainty band was not too much wrong since both 
effects were roughly compensating for each other. One way 
to quantify the discrepancy of Pb -I- Pb data from nor- 
mal nuclear absorption consists in fitting the points corre- 
sponding to most central collisions, i.e. beyond L = 8fm, 
with an exponential function of L (Fig. ||). The result is an 
effective additional absorption cross section of 9 mb, with 
uncertainties of 0.6 and 2.6 mb depending on whether one 
uses the "minimum bias" or the measured J/^ over Drell- 
Yan ratio. With respect to the reference absorption cross 
section of 5.8 ± 0.7 mb, the significance of this additional 
absorption amounts to 9.8 or 3.3 standard deviations, re- 
spectively, with a clear advantage to the "minimum bias" 
ratio, because it uses the normalization to the inclusive 
centrality distribution. 



6 Conclusion and perspectives 

A new data presentation has been proposed for results 
from nucleus-nucleus collisions. Its domain of application 
is not limited to J/W and Drell-Yan production processes 
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Fig. 5. "Minimum bias" J/<f' over Drell-Yan ratio divided by 
the normal absorption suppression as a function of the mean 
nuclear path length L (Fig. 9 of [H].) An uncertainty band 
for normal absorption, and an exponential fit for additional 
absorption in most central collisions beyond L = 8fm have 
been added to the original figure 



which have been chosen for illustration. For any "p" pro- 
cess one ends with its yield per nucleus-nucleus collision 
as a function of the estimated squared impact parameter. 
Since the normalization of the yield refers to the most 
probable, i.e. inclusive, process, one keeps the best statis- 
tical accuracy for each process. It seems to be a good way 
for going as far as possible with raw experimental data, 
sticking as closely as possible to them while trying to show 
physical quantities of interest. Could it be the best way to 
present experimental results concerning nucleus-nucleus 
collisions before comparison to any model? From the the- 
oretical side, one would like to compare experimental re- 
sults with results from model calculations on plots using 
the most relevant variable from the model, for instance the 
number of participants, the number of nucleon-nucleon 
collisions, the mean path length in nuclear matter, etc. 
The proposed data presentation could serve as the com- 
mon basis before going to any of these plots. 

For future nucleus-nucleus experiments that will take 
place at RHIC and LHC, the present work shows that it is 
possible to study any process independently of all others. 
The experimental results to be presented for each process 
have a direct physical meaning and are easily compared 
with model calculations. The only requirement is the mea- 
surement of the inclusive differential cross section with re- 
spect to at least one centrality variable used to sort out 
events according to impact parameter. It is essential that 
such inclusive measurements be available in experiments 
to be performed at RHIC and LHC. 
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